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Indirect search for Solar dark matter with AMANDA and IceCube 
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Abstract: It has been suggested that the Milky Way dark matter halo consists of Weakly Interactive Massive Particles 
(WIMPs) that scatter and accumulate in the center of gravitational wells like the Sun. In certain models this provides a 
source of high energy neutrinos originating in the centre of the Sun that may be detected by a neutrino detector on the 
Earth. We describe a search for Solar WIMP dark matter using AMANDA/IceCube data from 2008 with the 40-string 
IceCube detector configuration. We present combined limits using data acquired from 2001-2008 on the WIMP induced 
muon flux from the Sun for dark matter particle masses between 50 and 5000 GeV. We test WIMP model hypotheses of 
the lightest MSSM neutralino and the lightest Kaluza-Klein particle. 
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1 Introduction 

The intriguing question of why less than one fifth of the 
matter in our universe seems to be of the ordinary sort that 
emits electromagnetic radiation and the rest being invisible 
by all other means than gravitational inference, has been an 
outstanding problem in physics for almost 80 years. The 
most popular current hypotheses involve Weakly Interact- 
ing Massive Particles (WIMPs) that are distributed in and 
around galaxies as a halo. In several theories these particles 
will accumulate in massive objects, like the Sun, by scatter- 
ing (weakly) multiple times. With each scatter momentum 
is lost and the particle becomes gravitationally trapped in 
the centre regions of the object. Should the WIMPs then 
be Majorana particles they will pair-wise annihilate and, 
in some channels, give rise to a neutrino signal that is de- 
tectable on Earth. 

In this paper we present an explicit search for Solar WIMPs 
described by the Minimal Supersymmetric Standard Model 
(MSSM) [1], where the lightest supersymmetric particle 
xf, henceforth denoted x, is a promising dark matter can- 
didate. The two most extreme annihilation channels are 
studied; x ~^ bb producing the softest possible neutrino 
spectrum, and x W^W^ producing the hardest possi- 
ble spectrum. Note that for < mw we consider instead 
the channel x — ^ t~^t~ . 

The same event selection is also applied to search for an- 
other potential dark matter candidate, the Lightest Kaluza- 
Klein Particle (LKP), arising from theories of Universal 
Extra Dimensions (UED) [2]. The 5-dimensional model 
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Figure 1: IceCube, depicted from above, in its 40-string 
configuration (squares) together with AMANDA (circles). 
The dots show the full 86-string IceCube detector in 2010. 

we considered here is described by the mass m^^-* of the 
LKP, the first photon excitation, and the mass splitting 
A^(i) = (TOq(i) — TO^(i))/m^(i), where m^(i) is the mass 
of the first quark excitation. The Monte Carlo defining the 
signal here, with A^(i) = 0, was taken from [3]. 

The world's largest neutrino observatory, IceCube [4], was 
halfway to completion in 2008 with 40 strings, each con- 
taining 60 digital optical modules, deployed in the lay- 
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out shown in figure 1. AMANDA, the predecessor de- 
tector to IceCube, underwent a low-energy trigger upgrade 
and integration of its Data Acquisition System (DAQ) with 
that of IceCube. As seen in Figure 1, AMANDA, with a 
markedly lower energy threshold than IceCube, was hori- 
zontally enclosed by the 2008 IceCube array. Data from 
this year thus provided a unique opportunity to use the two 
detectors as one. In particular, IceCube could be used as 
an efficient veto for horizontal atmospheric muons enter- 
ing AMANDA. This opportunity was not missed, and we 
present here the result of a search for Solar dark matter us- 
ing this data, and follow up by combining the result with 
previous AMANDA and IceCube limits. 

2 Data and simulation 

Data 

An analysis was undertaken on the 2008 IceCube- 
AMANDA dataset amounting to a total Uvetime of 149 
days. The livetime considers the period when the sun was 
below the horizon (March to September). A triggered sam- 
ple of 1.7x10^" events exists in this dataset, the vast ma- 
jority of which are atmospheric muons. 

The chosen event-sample was to be reduced by a factor 
of ~10^^ using variable cuts and a machine learning al- 
gorithm in an effort to maximise the separation of signal 
and background. When optimising the event selection, de- 
scribed further in section 3, data from November 2008 was 
used and then discarded. 

Simulation 

The signal Monte Carlo was generated with DarkSUSY 
[5], simulating WIMPs and LKPs for a range of different 
particle masses. As mentioned in section 1, we generated 
two MC samples for each MSSM mass corresponding to 
soft and hard neutrino spectra. 

The atmospheric muon background was simulated with 
CORSIKA [6] using the Horandel [7] model for the cos- 
mic ray spectrum and composition. The background of at- 
mospheric was generated with the Honda et al. [8] con- 
ventional flux model and the Enberg et al.[9] prompt flux 
model. After simulating charged particle propagation with 
MMC [10] and Hght propagation in the ice with PHOTON- 
ICS [11], the detector response was finally simulated using 
the IceCube simulation software IceSim. 

3 Event selection 

This analysis divided the data into two data streams: One 
stream used IceCube as the main detector searching for 
high WIMP masses, using as a signal template the 1000 
GeV WIMP model with the hard spectrum. The second 
stream used AMANDA as the main fiducial volume and 




Figure 2: Effective volume for a 1000 GeV WIMP signal 
on the left y-axis and background rates on the right y-axis, 
as functions of the accumulated data reduction cuts in the 
IceCube-stream. In the steep step to the final level a cut 
was applied on SVM output values. 

IceCube as a sophisticated veto, with the 100 GeV soft- 
spectrum WIMP as a signal template. 

In each stream a series of linear cuts were applied to data 
and MC, optimised individually using the signal templates 
and data from November 2008 as the estimate for the back- 
ground. In the AMANDA-stream no more than 4 hits 
were allowed in IceCube, effectively favouring low-energy 
events. In both streams we made log-likelihood reconstruc- 
tions of the particle tracks, requiring them to be upward- 
going and good fits to the hit information. The progression 
of data reduction is visualised in figure 2. 

After applying the linear cuts a Support Vector Machine 
(SVM) [12] was trained to separate signal from back- 
ground. A separate signal MC sample and November data 
as background were again used, the both of which were dis- 
carded after the training. When applying a trained SVM on 
a set of events, a single number Q £ [0, 1] is obtained for 
each event that characterises it as signal-like or data-like, 
which can be used as a powerful multivariate cut parame- 
ter. 

Starting with a set of about 20 variables, an iterative 
method was used to eliminate variables that did not con- 
tribute significantly to the overall separation. In this way, 
12 and 10 variables were chosen for the SVMs in the 
AMANDA- and IceCube-stream respectively, covering a 
large part of the available phase-space and showing modest 
internal correlations. 

The final cut on Q was set at 0. 1 , which optimised the sensi- 
tivity of the dark matter search, i.e., the median 90% upper 
limit in a large set of random samples drawn from a back- 
ground distribution without any signal content. It is noted 
that the neutrino spectra from the Kaluza-Klein signal are 
very similar to the hard MSSM spectra once convoluted in 
the propagation through the Sun and processed as IceCube 
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triggers. For this reason, no separate optimisation was nec- 
essary to study this potential signal. 

Uncertainties 

A number of systematic effects were studied which impact 
the effective volume of the signal. Uncertainties in the 
neutrino oscillation parameters, neutrino-nucleon cross- 
section, muon propagation through ice and calibration of 
detector time and geometry, all change Vcff less than 5% 
when individually shifted la from their nominal values. 
The uncertainty in light propagation through the ice and the 
absolute sensitivity of the PMTs remain the dominant sys- 
tematic effect. This was estimated from dedicated Monte 
Carlo studies. The total systematic uncertainty, the squared 
sum of all effects, was calculated for 3 template energy re- 
gions: The AMANDA-region, the combined region and the 
IceCube-region were given the uncertainties 26%, 25% and 
15% respectively. 

4 Search for a Solar excess 

At the final cut level, the 2008 data sample contains a total 
of 3012 events. We examine the properties of the space an- 
gles (ipi) between the reconstructed track and the Sun. The 
angular resolution, measured as the median space angle for 
neutrinos in the signal MC sample, is 2.2° for the 1000 
GeV hard-spectrum signal and 9.0° for the 100 GeV soft- 
spectrum signal. The distribution of IceCube data in the 
direction of the Sun is shown in figure 3 together with the 
expected background. This background probability density 
function (PDF) is constructed from actual data using the 
reconstructed zenith-angles and random, oversampled, az- 
imuth angles. 

Since the AMANDA- and IceCube-streams have com- 
pletely non-overlapping event samples, we will treat them 
as separate sets of observations, denoted AOS and IC08, 
that can be combined using techniques described below. In 
addition to these samples, we also include a sample from 
an earUer analysis of AMANDA data from 2001-2006 [13] 
(AO 1-06, around 2000 events) and the sample from the 22- 
string IceCube data of 2007 [14] (IC07, 6946 events) to 
prepare a final analysis that incorporates up to four inde- 
pendent event samples. For all WIMP channels we com- 
pare the individual sensitivities and exclude, for compu- 
tational efficiency, samples yielding a sensitivity 10 times 
weaker than the strongest. In the results for the Kaluza- 
Klein models, e.g., AOS did not contribute significantly and 
was dropped in favour of IC07 and IC08. 

The two older samples were both obtained in a similar man- 
ner to what is described in this paper. In AO 1-06 a Boosted 
Decision Tree (BDT) was trained and optimised for each 
WIMP channel using 21 variables, in this way performing 
several parallel searches in a detector livetime of 812 days. 
The optimal BDT cut was chosen to minimise the number 
of events needed to have 90% probability of a 5ct discovery. 
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Figure 3: The distribution of space angles between the re- 
constructed event track and the Sun, for IceCube data at 
the final cut level and the expected background, in the re- 
gion near the Sun. The error bars on the background are Icr 
Poisson standard deviations. 

The IC07 sample resulted from a combination of straight 
cuts and two SVMs. Each SVM was trained with 6 differ- 
ent uncorrected variables, and the optimised cut was made 
on the product of the two output values. 

The goal here is to obtain a confidence interval on the num- 
ber of signal events, fj,, in the total data set. We use the 
same procedure and profile-likelihood method as outlined 
in [15], but adapt it for the more general case of TV indepen- 
dent experiments. The total likelihood of having /i signal 
events in the data set becomes 



N 



N 



(1) 



where each fij is weighted with the livetime and effective 
volume of the individual experiment using 



^livc^efF 



4-llVO 
fc \ 



(2) 



so that = ji. In equation 1, f{'4>\ii) is the total PDF 

of expected background plus a signal content of size fi. 

Table 1 shows the upper limit, at 90% confidence level, for 
fjL (or the entire confidence interval if the lower bound is 
not zero) together with the sensitivity as defined in section 
3, the best fit of fi to data and the probability of observ- 
ing at least the current deviation from the sensitivity. This 
was consistent with /i = 0, and the upper edge of the inter- 
val was converted to an upper limit, following [15], on the 
muon flux from the Sun. 
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29.2 
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23.7 
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1000 


14.0 
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17.1 
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0.21 


5000 
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17.1 
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0.20 
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50 


20.6 


16.2 


3.0 


0.33 


100 


19.4 


12.2 
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0.20 


250 


15.3 
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0.0 


0.43 


500 


10.9 


16.8 


0.0 


0.26 


1000 


6.8 


15.9 


0.0 


0.13 


3000 


9.3 


16.4 


0.0 


0.22 


5000 


8.7 


15.8 


0.0 


0.22 



Table 1 : With /i as the number of signal events in the final 
combined sample, the columns show the 90% CI or upper 
limit, the sensitivity, the best fit to data and the probability 
to observe at least this deviation from the sensitivity 



5 Results 

Taking into consideration the uncertainties discussed in 
section , figure 4 shows the 90% CL upper limits on the 
muon flux from the Sun by self-annihilating dark matter 
modelled by MSSM and Kaluza-Klein. In the MSSM plot 
we outline the MSSM model space not yet excluded by 
the direct detection experiments CDMS [16] and XENON 
[17], generated by an extensive parameter scan. The lower 
figure indicates the available Kaluza-Klein model space for 
considered values of A^(i) . We refer the reader to [3] for 
a more detailed account of probing the Kaluza-Klein space 
with IceCube. 

6 Discussion 

We have presented an indirect search for Solar WIMP 
dark matter using data taken with the integrated 
AMANDA/lceCube detector in 2008. Combining this data 
with two other event samples from AMANDA and IceCube 
from 2001-2007 we establish the strongest current upper 
limit on the muon flux from WIMP dark matter annihila- 
tion in the Sun. The completed IceCube detector has now 
started to take data and, together with the DeepCore low- 
energy extension, the expected sensitivity will probe sub- 
stantial regions of the previously untested WIMP parame- 
ter space. 
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1 Introduction 

While the presence of dark matter in the universe has been 
inferred through its gravitational interactions, it has yet 
to be directly or indirectly observed. One of the most 
promising and experimentally accessible candidates for 
dark matter are so-called Weakly Interacting Massive Parti- 
cles (WIMPs). In theories of the Minimal Supersymmetric 
Standard Model (MSSM), the WIMP can take the form of 
the lightest neutralino, x [1]^ while in the framework of uni- 
versal extra dimensions (UED) [2] it would be the lightest 
Kaluza-Klein (KK) particle (LKP). Here, we consider the 
first excitation of the KK photon, 7^^^, to be the LKP. In 
current models, WIMPs are predicted to have a mass in the 
range of a few lO's of GeV to a few TeV. Whatever their un- 
derlying physics, these WIMPs may be swept up by the Sun 
on its transit through the galactic halo and become grav- 
itationally botmd by scattering weakly on solar nucleons. 
Over time, this leads to an accumulation of dark matter in 
the center of the Sun, exceeding the mean galactic density. 
Self-annihilation to standard model particles may result in 
a flux of high energy neutrinos that will be spectrally de- 
pendent on the annihilation channel and WIMP mass, and 
which can be searched for as a point-Uke source with neu- 
trino telescopes such as IceCube. Several analyses [3, 4, 5] 
performed on the data from the IceCube detectors have al- 



ready set stringent limits on the WIMP induced muon flux 
from the Sun and the spin-dependent WIMP-proton scat- 
tering cross section comparable to or better than those of 
direct detection experiments. 

1.1 The IceCube Telescope 

The IceCube Neutrino Telescope [6] records Cherenkov 
light in the ice from relativistic charged particles created 
in neutrino interactions in the vicinity of the detector By 
recording the arrival times and intensities of these photons 
with optical sensors, the direction and energy of the muon 
and parent neutrino may be reconstructed. IceCube instru- 
ments 1 km^ of glacial ice at the South Pole with 5160 
Digital Optical Modules (DOMs) on 86 strings deployed 
between depths of 1450m-2450m. Eight more densely in- 
strumented strings optimized for low energies plus the 12 
adjacent standard strings at the center of the detector geom- 
etry make up the DeepCore subarray, which increases the 
sensitivity at low energies and substantially lowers the en- 
ergy threshold. In addition, by using the surrounding Ice- 
Cube strings as a veto, DeepCore will enable searches at 
low energies in the southem hemisphere, transforming Ice- 
Cube into a full sky observatory. The current analysis of 
data uses the 79 string configuration of the detector, includ- 
ing 6 DeepCore strings, marked with squares in Figure 1. 
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Figure 1: Top view of the final deployed IceCube-86 detec- 
tor configuration. Standard IceCube strings are shown in 
circles, initial DeepCore strings in squares, 'in-fill '-strings 
in DeepCore in triangles and their original design position 
with crosses. All strings without solid markers make up the 
79 string configuration. 

The sensitivity study we present for the full 86 string de- 
tector uses the initially proposed DeepCore geometry con- 
sisting of the same 6 additional strings, as present in the 79 
string detector. 

2 86 String Sensitivity Study 

In order to gauge the dark matter physics potential of Ice- 
Cube as well as the impact of the DeepCore subarray, we 
have performed a detailed study to determine the sensitiv- 
ity of the 86 string detector to signals originating from dark 
matter annihilations in the center of the Sun. (note, Deep- 
Core infill strings, as indicated in Fig. 1, are assumed at 
their original design position within this study.) Contrary to 
previous estimates, this study was performed as a full anal- 
ysis in all details, including detailed data processing and 
event selection, while making conservative choices where 
possible. This gives us a reaUstic expectation of the capa- 
bilities of IceCube to observe dark matter induced signals, 
given the state of data extraction, reconstruction, and signal 
discrimination techniques available at the time of the study, 
conducted prior to the start of 79 string data taking. 

2.1 Data Sample 

High statistics atmospheric muon and neutrino background 
samples are simulated including both single events as 
well as multiple muon coincidences within the detec- 
tor. The propagation of muons and photons in the ice is 
simulated [7, 8] taking measured ice properties into ac- 
count [9]. All signal simulations use DarkSUSY [10] and 
WimpSim [11]. 



The primary WIMP annihilation spectrum is extremely dif- 
ferent between x and 7^^\ For the LKP, branching ratios 
are fixed from theory, and therefore a 'true' spectrum can 
be simulated, whereas every MSSM model will result in a 
different set of branching ratios. To cover the full range, 
we assume 100% branching into the two end points of the 
spectrum, i.e the 'soft' bb and 'hard' W^W~ (t+t" be- 
low 80.2 GeV) channels. WIMP signals are generated for 
7 different WIMP masses with the 'hard', 'soft' and 'LKP' 
spectra. The selected masses for the Neutralino models 
range from 50 GeV to 5000 GeV and for the LKP models 
between 300 GeV to 1500 GeV. Three different algorithms 
are used to trigger the detector; the standard IceCube sim- 
ple majority trigger 'SMT' with multiplicity threshold 8, 
a string trigger (5 hit DOMs within 7 adjacent DOMs on 
the same string within 1000 ns) and a lower threshold SMT 
trigger for the fiducial DeepCore volume optinnized for low 
energies. 

2.2 Event Selection 

After performing likelihood-based reconstructions to 
achieve the best possible estimate of the muon direction, a 
series of increasingly stringent event selections are applied 
to remove the cosmic-ray induced backgrounds. Events are 
selected which appear to originate from the direction of the 
Sun and which exhibit low energy signatures. These cuts 
are not optimized for each model individually, but rather 
are chosen to provide good general signal retention, rang- 
ing from 25% to 70%, while rejecting as much background 
as possible (~ 4 x 10~^). The remaining signal and back- 
ground is separated through the use of multivariate train- 
ing algorithms. In order to determine the best possible 
approach, multiple algorithms were tested, including sup- 
port vector machines (SVM), neural networks (NN), and 
boosted decision trees (BDT). It was found that the best 
sensitivity could be achieved through the combined use of 
a BDT and a NN applied on independent sets of observ- 
ables with low correlation and high discrimination power. 
In order to further maximize the sensitivity across all mod- 
els, and taking into account the fact that observables are 
often strikingly different in different energy regimes, these 
machines were trained for 2 classes of models, 'high en- 
ergy' and 'low energy', depending on the resulting muon 
spectrum after dark matter annihilation and neutrino inter- 
action in the ice. A selection was made on the product of 
the outputs of the machines to maximize the Model Rejec- 
tion Factor (MRF) [12] of the analysis. This selection was 
then tightened in order to achieve a conservative result in 
light of the fact that the study is limited by statistics on the 
background simulation. With respect to the straight cuts 
level, the background reduction is 1 x 10~^ with signal ef- 
ficiencies ranging within 20% to 30%. 

From the two final 'low energy' and 'high energy' event 
selections of the signal simulation, we derive the effective 
area for muon neutrinos from the direction of the Sun as 
a function of neutrino energy, see Fig. 2. Also shown in 
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Figure 2: Effective area and median angular error for the 
final event selection as a function of neutrino energy in the 
range 20— 1000 GeV, for muon neutrinos from the direction 
of the Sun for the low and high energy cut selection. The 
result is an average over the austral winter. 

the figure is the median angular error between the recon- 
structed muon and the neutrino directions, 6. The result is 
an average over the austral winter, during which the Sun is 
below the horizon. 

2,3 Sensitivity 

After final event selection, the atmospheric muon back- 
ground reduction is < 5.2 • 10^^, which implies that sur- 
viving events in the background sample are dominated by 
atmospheric neutrinos. We search for an excess of neu- 
trino events over this expected background from the direc- 
tion of the Sun in a specifically determined search cone, 
the angular size of which is determined by maximizing the 
MRF= /if^Vn^, where Us is the number of simulated sig- 
nal events remaining in the search cone, and jfg^'^" is the av- 
erage Feldman-Cousins 90% confidence level upper limit 
on the signal given an expected background nbg and no 
true signal content. [13] We derive this 90% C.L. average 
upper limit, or sensitivity, for the 86 string detector con- 
figuration assuming a livetime of 180 days, and from this 
calculate the limit on the neutrino to muon conversion rate 
Fy^^. Using the signal simulation, we can convert this rate 
to a limit on the WIMP annihilation rate in the Sun, Fa. 

Assuming equilibrium between WIMP capture and anni- 
hilation in the Sun, the capture rate is then directly pro- 
portional to the WIMP-proton scattering cross sections a si 
and (TsD. By assuming that the capture is dominated either 
by spin-independent or spin-dependent scattering, it is pos- 
sible to derive sensitivities for a si and aso [14]. Results 
on <JSD are shown in Figure 3 for the case of the LKP and 
the neutralino, and compared to other experimental limits. 



3 79 String Data Analysis 

Data taking with the 79 string configuration of IceCube (in- 
cluding 6 DeepCore strings) took place between May of 
2010 and May of 2011, and analysis of this dataset is un- 
derway. Several improvements have been realized since 
the completion of the study detailed above and we describe 
their potential impact on the sensitivity of the analysis here. 
While this work is ongoing, it is anticipated that results will 
be seen in the coming months. 

• The low energy trigger in actual data taking was reduced 
from the planned 4 required hit DOMs within the Deep- 
Core fiducial volume to a less stringent 3 hit DOMs. This 
further extends the acceptance of the lowest energy events, 
and boosts the sensitivity for low mass WIMPs. 

• In previous datasets, hit information was recorded only 
for those DOMs which had a timing coincidence with their 
vertical neighbors. This virtually eliminates random noise 
hits, but at the same time removes a not insubstantial frac- 
tion of real physics hits. Beginning in 2009, this condition 
is relaxed, and limited hit information is read out for all 
DOMs which see light in a triggered event. Noise clean- 
ing algorithms have been developed to clean this expanded 
dataset of unwanted noise while retaining physics informa- 
tion vital for low energy analysis. 

• Rather than considering only those events in a search 
cone around the angular position of the Sun, it is possible 
to conduct the search utilizing a more sophisticated likeli- 
hood technique. Probability distribution functions (PDFs) 
for WIMP signal and Vatm background can be constructed 
relative to the solar position. The signal PDF, is well 
known for each model from simulation, while the back- 
ground PDF, fbg{^), can be found by randomizing the az- 
imuth angle in the final event sample. The likelihood that 
an event belongs to either of the PDFs can be evaluated, and 
the most likely signal contribution in the final event sample 
determined. This method is superior to a search cone in 
that all events are considered, while stiU penaUzing those 
that lie further away from the Sun's location. 

• Extending searches to the southern hemisphere opens up 
the entire period when the Sun would otherwise be above 
the horizon, or the austral summer. This doubles the po- 
tential Uvetime of the experiment, though the challenges 
in reducing the downgoing atmopsheric muon backgrounds 
are immense. Much of this background can be removed by 
considering DeepCore and its surrounding strings as a fidu- 
cial volume and using the surrounding IceCube strings as 
an active veto. Furthermore, techniques have been devel- 
oped to determine the UkeUhood that a reconstructed muon 
originates within the detector volume, a sign that it was the 
result of a neutrino interacting in the ice. 

Taking these and other improvements into account, we ex- 
pect to be able to substantially improve the sensitivity of 
the analysis over the basehne described in the study above, 
most notably at the lowest energies, where background dis- 
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Figure 3: (left plot) Limits on the LKP-proton SD scattering cross-section from IceCube [3] (dashed line) and the final 
sensitivity for IceCube-86 (solid line) compared to the theoretically allowed region of m^(i) and A^ii) . The region below 
TO^(i) = 300 GeV is excluded by collider experiments [15]. The upper bound on m^(i) corresponds to the over-closure 
limit for each individual LKP model [16]. The darker regions (blue in on-line version) indicate the overlap with two 
different ilcDM^i^ intervals [17]. (right plot) Sensitivity to x on proton spin-dependent cross section at 90% confidence 
level with respect to x mass. IceCube current best limits [3] are shown in solid markers, IceCube-86 sensitivity is shown 
in non-filled markers and compared to other experiments. Soft WIMP models are indicated by the dashed lines, whereas 
hard models are shown in solid lines. The indicated model space represents a scan over the allowed MSSM parameter 
space, accounting for all current experimental constraints. 



crimination is difficult and well reconstructed signal event 
rates are low. 

4 Outlook 

Construction of the IceCube telescope was completed in 
December of 2010, creating a sensitive platform for the 
discovery of Dark Matter candidates in the energy range 
of 50 GeV to 1 TeV. The co-located DeepCore subarray 
pushes the energy threshold of IceCube into the 10 GeV 
regime and allows for the possibility to extend searches to 
the southern hemisphere. For searches of dark matter an- 
nihilation in the center of the Sun, this effectively doubles 
the experimental livetime while greatly improving the sen- 
sitivity for low mass WIMP models where the bulk of the 
neutrinos and daughter muons are at these very low ener- 
gies. 

We have conducted a detailed simulation-based sensitivity 
study utilizing the full 86 string detector geometry to better 
understand the possible physics reach of the detector and 
have shown that with only 180 days of livetime, we will 
improve the limits on the spin-dependent WIMP-nucleon 
scattering cross section by an order of magnitude. In addi- 
tion, this study shows that the inclusion of the DeepCore 
detector significantly improves the reach at low WIMP 
mass, allowing IceCube to probe the full range allowed by 
current models. 

Data-taking with the 79 string configuration of the detector 
has concluded in May of 2011 and analysis of this dataset 
has begun. Many improvements in the detector configura- 
tion, extraction of useful information from the data, and in 



analysis techniques for maximizing the WIMP signal lead 
us to believe that we can substantially improve upon the 
predictions described in the above sensitivity study. With 
these improvements in mind, in the coming months we will 
be able to set the most stringent limits to date on the WIMP- 
proton spin-dependent scattering cross section, rule out in- 
dividual models, or potentially detect an excess neutrino 
flux from the Sun, strongly indicating the capture and an- 
nihilation of dark matter in its core. 
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1 Introduction 

The existence of a cold, non-baryonic dark matter compo- 
nent of the Universe is implied by a variety of astronomi- 
cal observations. The velocity distribution of galaxy clus- 
ters suggests a mass content significantly higher than what 
should be expected from observation of the luminous mass. 
On smaller scales, the rotation profiles of galaxies show 
discrepancies, which hint at the existence of a dark matter 
halo reaching out far beyond the baryonic disc and bulge. 
Based on e.g. the anisotropy in the cosmic microwave back- 
ground observed by WMAP [1], a determination of cos- 
mological parameters is possible, which yields the concor- 
dance model of cosmology, with a dark matter content of 
about 22%. Large scale structure formation and N-body 
simulation of cosmic evolution require the existence of cold 
dark matter, which is expected to consist of massive non- 
relativistic particles. Weakly Interacting Massive Particles 
(WIMPs) form a class of promising dark matter candidates, 
that are expected to have masses from a few GeV up to 
several lOOTeV [2, 3]. Assumed to be relics from the Big 
Bang they would naturally yield correct cosmic relic abun- 
dances. A promising WIMP candidate is the neutralino, the 
Ughtest stable particle introduced by the hypothetical min- 
imal super-symmetric extension to the Standard Model. 

WIMPs, if assumed to be Majorana particles, can self- 
annihilate and produce a variety of Standard Model par- 
ticles in the final state. The dark matter annihilation in 
the Milky Way would hence yield an excess flux of neutral 
messenger particles (7, i^) from the direction of the Galac- 



tic Center, where the dark matter density is expected to 
peak. Further, it would also lead to a large-scale anisotropy 
over the sky. Such a flux of annihilation products can be 
detected by earth-based, air-borne or satellite experiments, 
or the absence of such a flux can be interpreted in terms 
of upper Umits on the velocity-averaged self-annihilation 
cross-section, (itav) , of dark matter particles in the con- 
sidered model. In this paper we present two approaches to 
the search for dark matter with IceCube; an analysis of a 
possible large-scale anisotropy and a search for an excess 
flux from the direction of the Galactic Center. 



2 Neutrino Flux from Dark Matter Annihi- 
lation in the Galactic Halo 

Dark matter density profiles based on observations of dark 
matter dominated galaxies tend to have a rather flat density 
distribution in the central region, while profiles based on 
N-body simulations show a steep, often divergent increase 
of density towards the Galactic Center [4, 5]. An often- 
used parametrization of (sphericaUy symmetric) dark mat- 
ter profiles is 



Po 



"\ (/3-7)/« 



(1) 



where, r is the distance from the Galactic Center and 
is the scaUng radius. The commonly used NFW-profile is 
then defined by the parameters: r, = 20kpc, (a. /I 7) = 
(1,3, 1) [6]. Po is chosen so that the locally (at radius 
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Figure 1 : The NFW, Kravtsov, and Moore halo profiles are 
compared. 




Figure 2: J(^) is shown for the NFW, Moore, and 
Kravtsov profile. 

of the solar circle i?sc) assumed dark matter density is 
0.3 GeV/c^cm~^. Figure 1 illustrates the NFW profile 
along with two other profiles with different scaling param- 
eters; the flat Kravtsov profile [7] and the cuspy Moore pro- 
file [8]. 

Following [9], the expected neutrino flux from WIMP an- 
nihilation depends on the squared dark matter density inte- 
grated along the line of sight 



J -RscPsc 



where 'J is the opening angle between the line of sight and 
the Galactic Center and i?sc and p^^ are used as scaling pa- 
rameters, so that Ja(^') is dimensionless. The upper line- 
of-sight integration limit , Imax, is chosen to be approxi- 
mately the Milky Way radius. Figure 2 shows the shape of 
Ja(vE') for the above mentioned profiles. The differential 
flux can be written as: 



1e 



2 



RacPsc dM 



47rm^ 



dE' 



(3) 



The factor I/Att comes from isotropic emission, is the 
mass of the dark matter particle, and ^ is the energy spec- 
trum of the final state neutrinos. For the presented analyzes, 
DarkSUSY [10] scans have been performed, for a variety 
of benchmark WIMP masses and annihilation channels. 



3 IceCube Neutrino Observatory 

IceCube is a km'^-size neutrino detector located at the ge- 
ographic South Pole. Construction has been completed in 
December 2010. It consists of 86 strings instrumented with 
60 Digital Optical Modules (DOMs) each, including a low 
energy sub-array DeepCore, that has an energy threshold 
of the order of lOGeV. The DOMs are located at a depth 
ranging from 1.45 km to 2.45 km in the Antarctic glacier. 
The IceCube strings are arranged in a hexagonal symmetry, 
with an inter-string spacing of 125 m and a DOM distance 
of 17 m. Eight densely spaced DeepCore strings are de- 
ployed between regular IceCube strings and have a DOM 
spacing of 7 m and an inter-string spacing of 72 m. For 
the presented analyzes the 22- and 40-string configura- 
tions were used, which did not have any DeepCore strings. 

4 Galactic Halo Analysis with IceCube-22 

Between June 2007 and March 2008, the partiaUy- 
deployed IceCube detector was operated in a configura- 
tion with 22 strings and acquired 275 days of live-time, 
that were searched for neutrino signals from WIMP anni- 
hilations in the outer Galactic halo (Northern hemisphere). 
The sample contained 5114 up-going muon neutrino can- 
didate events, covering —5° and 85° in reconstructed decli- 
nation, with a purity of about 90% [11]. Despite the Galac- 
tic Center (RA 17h45m40.04s, Dec -29°00'28.1") being 
located on the southern hemisphere and not accessible in 
this sample, we expect an anisotropy to be present on the 
northern hemisphere. The anisotropy is caused by the dark 
matter distribution in the Milky Way halo, which would 
result in a larger neutrino flux from the region closer to 
the Galactic Center compared to that further away from it. 
We have searched for a neutrino anisotropy by comparing 
fluxes from an on-source region centered around the same 
RA as the Galactic Center to that from an equally sized 
off-source shifted by 180° (see figure 3). The analysis was 
performed in an unbiased fashion, by only counting events 
in the regions after their size and the analysis procedure 
were optimized using simulations only. 

We observed compatible numbers of events in the on- and 
off-source region: 1367 and 1389, respectively [12]. As we 
perform a relative comparison on data the uncertainties on 
the background estimate can be kept small. A systematic 
uncertainty of 0.3% due to any pre-existing anisotropy in 
the data caused by exposure and cosmic-ray anisotropy re- 
mains. The signal acceptance uncertainty is approximately 
30% and is dominated by the incomplete modeling of ice 
properties and limitations in the detector simulation. Us- 
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Figure 3: Relative expected neutrino flux in the northern 
hemisphere from self-annihilation in the Milky Way halo. 
The on-source region (solid line) is centered around largest 
the flux expectation at ARA = 0, while the off-source 
region is shifted by 180° in RA. 

ing equations (2) and (3), a limit on the self-annihilation 
cross section has been calculated and is shown in figure 4 
compared with the limits from the Galactic Center analy- 
sis, described in the next section. As the analysis uses the 
outer halo, the uncertainty on the choice of halo model is 
small as indicated by the error band on the limits. 

5 Galactic Center Analysis with IceCube-40 

The 40-string configuration of IceCube was taking data 
from April 2008 to May 2009, yielding a total detector live- 
time of 367 days. 

The highest neutrino flux from WIMP annihilation is ex- 
pected to come from a relatively wide region centered at 
the direction of the Galactic Center which, at the location of 
IceCube, is always about 30° above the horizon. Data from 
this direction is dominated by atmospheric muons, there- 
fore this analysis is based on the identification of events 
with an interaction vertex inside the detector (atmospheric 
muons produce incoming tracks) and it relies on the on- 
source/off-source method; based on Monte Carlo simula- 
tions, the width of a declination band (centered at the lo- 
cation of the Galactic Center) is optimized to maximize 
signal/ ^/background, assuming the NFW-profile. In this 
declination band, a window in right ascension is optimized. 
The optimum window sizes both in right ascension and 
declination were found to be ±8°. After correction for 
uneven exposure, as well as signal quality cuts, the un- 
certainty on the background prediction is reduced to the 
0.1%-level. Based on the above mentioned background es- 
timation, the expected number of background events in the 
signal region was 798819. The number of observed events 
was 798842. The difference of 23 events is compatible with 
the null-hypothesis, therefore a 90%C.L. -limit on the num- 
ber of signal events has been calculated (1168), following 
the Feldman-Cousins approach [13]. Using equations (2) 
and (3), a limit on the self-annihilation cross-section has 
been calculated and is shown in figure 4 along with the lim- 




[GeV] 

Figure 4: 90% C.L.-limits on the (cr^iv) from the IceCube- 
22 halo analysis (blue-shaded lines) [12], and the limits ob- 
tained from the IceCube^O Galactic Center analysis (sim- 
ple lines). For both analyses the lines from top to bottom 
correspond to the bb, W+W^, n^fJ.^ and i/D annihilation 
channels. The IceCube-40 limits are preliminary. 

its from the previous analysis. Figure 5 shows the obtained 
limits for the r channel, compared to the PAMELA/Fermi 
regions [14]. 

The IceCube^O limits are preliminary, since they do not 
include signal acceptance systematic uncertainty due to op- 
tical ice properties. 

6 Outlook on the Galactic Center Analysis 
with IceCube-79 

For IceCube-79, a dedicated Galactic Center data filter has 
been implemented and was taking data from June 2010 to 
May 2011. The filter consists of two parts. A so-called 
high energy part accepts all events with a reconstructed ar- 
rival direction within an angular window of ±10° in decli- 
nation and ±40° in RA with respect to the direction of the 
Galactic Center and if their brightness exceeds a zenith- 
dependent threshold. The so-called low-energy part ac- 
cepts events from a 15° wide zenith band around the Galac- 
tic Center, but applies a pre-scale factor of 3 on events from 
the zenith band, which have a distance of more than 20° to 
the Galactic Center in right ascension. Further restrictions 
for the low energy filter are a top veto defined by the upper 
5 DOMs, in which no hits are allowed, and a side veto. The 
side veto consists of the outer layer of IceCube strings; the 
earliest pulse is not allowed in this veto region. These fil- 
ter conditions allow for a preselection of tracks, which ap- 
pear to start within IceCube. Figure 6 shows a comparison 
of the effective area at filter level for IceCube^O and for 
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Figure 5: IceCube 90% C.L. upper limits on the ((Tav) 
from the Galactic halo with the 22-string and the Galactic 
Center with the 40-string array compared to the preferred 
regions for PAMELA data, and the region including Fermi 
data for annihilation to r+r^ [14]. 

the low and high-energy filter of IceCube-79. Especially 
in the low energy region below 100 GeV, more events are 
accepted. This improvement can be attributed to the Deep- 
Core array. If these events can be retained throughout the 
analysis cuts, considerable improvement is to be expected 
for exclusion limits on the self-annihilation cross-section in 
the low energy region. 

7 Conclusion 

Data collected with the partially instrumented IceCube 
neutrino detector has been searched for dark matter self- 
annihilation signals. Two independent analyses, target- 
ing the Galactic halo and Galactic Center, have been per- 
formed and resulted in observations consistent with back- 
ground expectations. Based on these results the dark mat- 
ter self-annihilation cross section was constrained to ^ 
lO^^^cm^s^^ for WIMP masses between 200 GeV and 
10 TeV for annihilation into T+r~ and /i^/i^. For a neu- 
trino line spectrum XX ^ annihilation cross sections 
larger than ^ can be excluded, assuming the 

NFW-profile for the Galactic Center analysis. Limits from 
the halo analysis are less halo-profile dependent, since the 
different models show similar behavior for larger distances 
from the Galactic Center. Despite the small dataset and 
less than half of the full IceCube detector, the limits al- 
ready probe a region of interest. A new dedicated filter 
stream for neutrinos from the Galactic Center has been im- 
plemented, that led to an increase in neutrino effective area 
at filter level of about two orders of magnitude at energies 
below 100 GeV. With the IceCube detector completed and a 
dataset available that is already more than three times larger 




Figure 6: Effective area for IceCube^O, and the two parts 
of the Galactic Center filter for IceCube-79 at online filter 
level. 

than the ones used for the presented analyzes, we expect 
to probe dark matter self-annihilation cross sections below 
^ Further, the Galactic halo analysis is cur- 

rently pursued using the DeepCore detector and the cas- 
cade channel (z^e, Vt)- It utilizes the excellent atmospheric 
muon veto capabilities with IceCube/DeepCore and lower 
atmospheric neutrino background in this channel. As the 
analysis targets a large scale anisotropy, the poor angular 
resolution of cascade events does not effect this analysis in 
a strong manner, and will allow for a further improvement 
in sensitivity. 
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1 Introduction 

The observational evidence for the existence of dark matter, 
from galactic to cosmological scales, is strong. However, 
its underlying nature remains unknown. A variety of the- 
ories provide candidate particles for cold dark matter [1]. 
Supersymmetry and Extra Dimensions predict new parti- 
cles with masses around the electro-weak scale, including 
a stable or long lived weakly interacting massive particle 
(WIMP) [1,2]. Such WIMPs form ideal dark matter candi- 
dates, whose masses are predicted to be in a range between 
a few tens of GeV to several TeV. 

A variety of standard model messenger particles (including 
high energy neutrinos) are expected to be produced as a re- 
sult of the self-annihilation or decay of WIMPs. The neu- 
trinos can be detected by high energy neutrino telescopes, 
while gamma-rays can be probed with ground based (e.g. 
H.E.S.S [3]) and space based (e.g. FERMI [4]) telescopes. 
Electron positron fluxes are sensitive to dark matter annihi- 
lations in the vicinity of the Sun (< 1 kpc) and have so far 
resulted in inconclusive results (e.g. PAMELA [5]). 

The dark matter self-annihilation cross section can be 
probed by looking for neutrino signals from the Galactic 
Center, Galactic halo [6] or dwarf galaxies. 

In this paper we discuss the first IceCube search for neu- 
trino signals produced by annihilating dark matter in dwarf 

galaxies surrounding the Milky Way. These objects are 
attractive targets for such searches because of their close 



proximity, relatively compact nature, and their large frac- 
tion of dark matter. The search wiU be used to probe the 
self-annihilation cross section by constraining the product 
of cross section and velocity averaged over the dark matter 
velocity distribution, {(jv), to probe the lifetime, r, and to 
make comparison with recent measurements performed by 
the FERMI collaboration. Sensitivities will be given for a 
set of selected benchmark annihilation chaimels for indi- 
vidual dwarf galaxies. 

2 Neutrino detection with IceCube 

The IceCube Neutrino Observatory [7], located at the geo- 
graphic South Pole, consists of the IceCube neutrino tele- 
scope and the IceTop air shower array. In the ice, a volume 
of one cubic kilometer of antarctic ice is instrumented with 
5160 digital optical modules (DOMs) deployed at depths 
between 1450 m and 2450 m. The DOMs are distributed 
over 86 electrical cable bundles that handle power trans- 
mission and communication with electronics located on 
the surface. Each DOM consists of a 25 cm Hamamatsu 
R7081-02 photomultiplier tube connected to a waveform 
recording data acquisition circuit. 

IceCube is sensitive to all flavors of neutrinos through 
Cherenkov light emission from secondary particles created 
when a neutrino interacts in the ice. Muon neutrinos are of 
particular interest since their extended track-like signature 
makes them relatively simple to identify and to reconstruct 
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Source 


Right asc. 


Declination 


Distance (kpc) 


Mass (10 Mq) 


J{lQ''^GeV'^/cm}') 


Segue 1 


10 07 04 


+16 04 55 


25 


1.58 


1 oc+3.Yt) 
-'-•^D_0.94 


Ursa Major II 


08 51 30.0 


+63 07 48 


32 


1.09 


n, ro+0.91 
U.OO_o.35 


Wilknan 1 


10 49 22.3 


+51 03 04 


38 


0.77 




Coma Berenices 


12 26 59 


+23 55 09 


44 


0.72 


16+°-^^ 

^•-■-0-0.08 


Ursa Minor 


15 09 08.5 


+67 13 21 


66 


1.79 


"•"*-0.18 


Draco 


17 20 12.4 


+57 54 55 


80 


1.87 


1 9+0.31 


Ursa Major I 


10 34 52.8 


+51 55 12 


106 


1.1 


Hercules 


16 31 02 


+12 47 30 


138 


0.72 




Canes Venatici II 


12 57 10 


+34 19 15 


151 


0.7 




Leo II 


11 13 29.2 


+22 09 17 


205 


1.43 




Canes Venatici I 


13 28 03.5 


+33 33 21 


224 


1.4 




Leo I 


10 08 27.4 


+12 18 27 


250 


1.45 




Leo T 


09 34 53.4 


+17 03 50 


417 


1.3 





Table 1: List of dwarf spheroidal galaxies ordered by distance from the Earth [14]. J values from [8] and for Segue 1 from 
[15]. The integral in Eq. 2 is performed over a circle of 0.25° radius for Segue 1 and 0.5° for the other sources. 



their direction with a few degrees precision at the detection 
threshold of 50 GeV. This analysis is based on data taken 
from May 2009 to May 2010 with an intermediate con- 
struction stage of the in-ice detector with 3540 DOMs on 
59 strings monitoring ~ 2/3 of a Gton of antarctic ice. The 
trigger rate in this configuration is between ^ 1600 Hz and 
~ 1900 Hz and reflects the seasonal modulation of muons 
produced in cosmic air showers. In future analyses, the 
denser low energy extension DeepCore, finalized in 2010, 
will considerably improve the detection of muons with en- 
ergies below a few hundred GeV. It will also provide sen- 
sitivity to celestial objects in the South by employing the 
veto capacity of the surrounding IceCube strings. 

The primary background in the search for neutrinos origi- 
nates from cosmic ray hadronic air showers produced in the 
Earth's upper atmosphere. The decay of pions and kaons 
results in a continuous stream of neutrinos and muons. 
High energy muons are capable of travelling long distances 
through matter before they eventuaUy decay, resulting in a 
down-going muon flux at the IceCube detector. We there- 
fore restrict ourselves to dwarf galaxies in the northern 
hemisphere, using the Earth as a shield. 

This paper discusses a sensitivity study of WIMP annihi- 
lation in dwarf galaxies. One of the appeals of this analy- 
sis is the ability to perform a direct comparison with the 
search for photons originating in a similar set of dwarf 
galaxies. We therefore follow the procedures outUned in 
a recent publication by the FERMI collaboration which set 
flux limits on several dwarf galaxies [8] and presented a 
prelinninary stacking analysis [9]. 

3 Dwarf galaxies target selection 

The number of identified dwarf galaxies has risen in re- 
cent years due to a systematic search by the SLOAN Dig- 
ital Sky Survey (SDSS) [10, 11]. Table 1 Usts galaxies in 
the Northern hemisphere that are accessible to IceCube in 



the up-going neutrino event sample. The expected neutrino 
flux from annihilating dark matter is [12] 



d<^{/\n,E) (av) dN 



dE 



87rm| dE 



(1) 



where denotes the WIMP mass and ^ the energy 
spectrum of the produced neutrinos per annihilation. The 
"J factor" is the Une-of-sight integral of the squared dark 
matter density: 

J{Afl) = [ dn ( p\{s)ds. (2) 

J AO. Jl.o.s. 

Note that this factor is highly dependent on the assumed 
dark matter distribution as the annihilation is propor- 
tional to the square of the dark matter density p^. For 
this paper we assume the commonly used Navarro-Frenk- 
White (NEW) [13] profile to describe the dark matter den- 
sity distribution. 

4 Analysis procedure 

In order to cover a large range of possible WIMP signa- 
tures, we study WIMP masses in the 0.3-10.0 TeV range 
and assume a few selected benchmark annihilation chan- 
nels: XX ~^ M"""/^"' T+r~, W+W~, bb, and lyp. Neu- 
trinos will have undergone extensive mixing through vac- 
uum oscillations over the distances travelled from the dwarf 
galaxies to the Earth. We determine neutrino flavor oscil- 
lations in the long baseline limit. The Monte Carlo results 
and selection efficiencies are obtained by reweighting the 
Monte Carlo sample with the expected muon energy dis- 
tributions for the WIMP masses and annihilation channels 
studied. 
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Figure 1 : Optimization of selection parameter for Segue 1 
and 1 TeV WIMPs annihilating to r+r^. Upper plot: Sig- 
nal mean upper limit (//go) divided by signal efficiency 
(Csig) versus bin size at an optimal cut value of 0.26. Lower 
plot: Signal mean upper limit divided by signal efficiency 
versus cut value at an optimal bin size of 1.8°. 

5 Event selection and sensitivity determina- 
tion for individual sources 

A pre-selection at the South Pole for up-going recon- 
structed muon tracks and two further likelihood reconstruc- 
tion and filtering steps reduces the data rate to ^ 2 Hz. We 
select events with reconstructed arrival direction 9 degrees 
below the horizon, as the lowest dwarf galaxy in the North 
is Leo I at a declination of 12°18'27". At this selection cut 
level the data samples are still dominated by atmospheric 
muons that are misreconstructed. The sample is then sepa- 
rated into one branch for each WIMP mass and annihilation 
channel. The event selection is optimized for the expected 
muon neutrino fluxes obtained with DarkSUSY [16]. Us- 
ing the Multivariate Analysis package of ROOT (TMVA) 
several Boosted Decision Trees (BDTs) are trained on sim- 
ulated signal and background samples, where the signal is 
weighted to the according energy spectrum for each Tree. 
The input variables include the likelihood of the track re- 
construction, the angular difference between different re- 
construction methods, an estimator for the individual an- 
gular resolution and the likelihood difference to a recon- 
struction that is forced to be down-going. 

Given the ^ 1 degree angular resolution of IceCube, the 
dwarf galaxies can be considered point-like. For each 
WIMP annihilation channel and dwarf galaxy, the event 
variable cuts and the circular area around the source loca- 
tion are optimized for the most stringent mean upper limit 
assuming no signal divided by the signal efficiency with re- 
spect to the pre-selection. Figure 1 shows the optimization 
for WIMPs of 1 TeV and assuming self-annihilations into 
T+T^. The background is estimated from a zenith band 
of ±2.5° around the source position. A Feldman-Cousins 
confidence interval construction [17] is used for the deter- 
mination of the upper limit. From the corresponding flux 



Figure 2: Expected mean upper limit on the neutrino flux 
versus WIMP mass for different dwarf spheroidals. Values 
for the T+T^ annihilation channel are shown. 
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Figure 3: Expected mean upper limit on the neutrino flux 
versus WIMP mass for different annihilation channels. Val- 
ues for Segue 1 are shown. 

a sensitivity for the WIMP annihilation cross section can 
be derived by using Eq. 1 . The average upper limit on the 
cross section is shown in figure 4. 

Figure 2 and 3 show the expected average upper limit on the 
neutrino flux as function of WIMP mass for different dwarf 
galaxies and annihilation channels, respectivly. The limits 
are based on an integrated lifetime of 334.5 days of data- 
taking with the 59-string detector configuration of IceCube. 

6 Source stacking 

The sensitivity can be improved by analyzing multiple 
sources simultaneously, a method called source stacking. A 
prerequisite of this technique is that the flux from all targets 
differs only in normalization. This is clearly the case for 
the WIMP annihilation signal of dwarf galaxies. Adding 
more sources to the analysis lets the signal grow linearly 
while the background grows proportionally to ^/N . Fig- 
ure 4 shows how the sensitivity on the annihilation cross 
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Figure 4: IceCube sensitivities for different number of 
stacked sources for one year of data in the r+r^ channel 
as function of the WIMP mass. 
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Figure 5: IceCube sensitivity (expected mean upper limit) 
for one year of data in the r+r^ channel as function of 
the WIMP mass. The sensitivity of the stacking analysis 
(Draco, Segue, Ursa Minor and Ursa Major II) is compared 
to results from MAGIC [18] and a preliminary combined 
limit from Fermi observations of 10 Dwarf Galaxies [9] 

section in the r+r^ channel improves while increasing the 
number of stacked sources. In Figure 5 we compare the 
IceCube sensitivity for four stacked sources assuming dark 
matter annihilation into r+r^, compared to results from 
MAGIC [18] and a preliminary combined limit from Fermi 
observations of 10 dwarf galaxies [9]. 

7 Conclusion 



data of these sources is collected continuously and extend 
to higher WIMP masses. The now operational full IceCube 
detector will be able to improve the sensitivity to dwarf 
speroidals further. 
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We have presented sensitivities for IceCube in the 59- 
string configuration for the observation of neutrino sig- 
nals from self-annihilating dark matter in dwarf spheroidal 
galaxies on the northern hemisphere. Sensitivities for the 
dark matter self-annihilation cross section are better than 
10~^°cm^s^^, for WIMP masses in a range of 300 GeV to 
several Te V. Searches are complementary to 7-ray observa- 
tion in most channels, but have the advantage that IceCube 
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1 Introduction 

The existence of magnetic monopoles was first suggested 
by Pierre Curie in 1894 [1]. However, a firm theoreti- 
cal grounding was not laid down until 1931 when Paul 
Dirac [2] demonstrated that magnetic monopoles are con- 
sistent with quantum mechanics. Dirac found that mag- 
netic charge is given by g = Ne/2a ~ 68. 5e, where a 
is the fine structure constant, and e the fundamental elec- 
tric charge. In 1974 t'Hooft and Polyakov independently 
described magnetic monopoles that are regular solutions of 
the field equations in certain groups of Grand Unified The- 
ories (GUT) and match the charge of the Dirac monopole 
[3, 4]. Within GUT the masses of magnetic monopoles can 
be estimated to be M oc A/a, where A is the unification 
energy scale of the theory. This results in a mass range from 
10^ GeV to 10" GeV for various GUT models. Because 
of these large masses magnetic monopoles are generally 
assumed to be relics of the early universe where they have 
been produced via the Kibble mechanism [5]. 

Analogous to electric charges, which are accelerated along 
electric field lines, magnetic monopoles are accelerated 
along magnetic field lines. The kinetic energy gained by 
a monopole is Ek oc gB^, where B is the magnetic field 
strength, and ^ is the coherence length of the field [6]. Dur- 
ing the Ufetime of the universe, reUc monopoles should 
have encountered enough accelerators to reach kinetic en- 
ergies of ^ 10^"' GeV. Thus monopoles with masses less 
than ~ 10^'' GeV should be relativistic. IceCube is able 
to detect magnetic monopoles travelhng through the de- 



tector at velocities greater than the Cherenkov threshold 
i(3 > 0.76). The radiation emitted by the monopole is pro- 
portional to (gn)^, where n is the index of refraction of the 
ambient medium [7]. Thus, in ice (n w 1.3) a monopole 
will emit ~ 8000 times more light than a bare muon of the 
same velocity. 

2 Detector 

The IceCube detector is a cubic kilometer-scale neutrino 
telescope. In its now completed state, IceCube consists of 
86 strings of 60 Digital Optical Modules (DOMs), each 
spaced out in a hexagonal pattern and deployed between 
1450 to 2450 meters below the Antarctic ice surface. For 
the data presented, we use the configurations of IceCube as 
of 2007 and 2008 with 22 and 40 strings respectively. Each 
DOM is configured to detect and digitize photon signals 
via a Photomultiplier Tube (PMT) and two waveform dig- 
itizers, a fast Analog to Digital Converter (fADC) and an 
Analog Transient Waveform Digitizer (ATWD) [8]. The 
ATWD has a sampling rate of ~ 300 Megasamples/second 
with a total of 128 samples and digitizes the incoming 
waveform across 3 channels representing different gain val- 
ues. The fADC runs with a sampUng rate of 40 Megasam- 
ples/second and can read up to 256 samples. This con- 
figuration allows the digitalization of short signals with a 
high dynamic range as well as long signals, although with 
a smaller dynamic range. 
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A DOM is triggered when the PMT signal exceeds a cer- 
tain threshold. However, the data is only read out in the 
case of a local coincidence with one neighboring DOM on 
the same string in a time window of 1/xs. This serves to 
reduce background originating from PMT noise. Since the 
waveforms produced by a monopole are expected to have 
a long time scale, the fADC provides greaterdistinction be- 
tween signal and background. Hence, the search strategies 
in this study mainly rely on the data provided by the fADC. 

3 Signal and Background Simulation 

The simulation of relativistic magnetic monopoles is done 
in three stages. Magnetic monopoles are generated uni- 
formly on a disk with their direction perpendicular to the 
disk. The disk itself is located ~ 1 km from the center 
of the detector pointing towards it at randomized orienta- 
tions. The radius of the disk is set to 650 m for the 22 
string configuration and 850 m for the 40 string configu- 
ration. Datasets were generated for four different speeds, 
13 = 0.995, /3 = 0.9, (3 = 0.8, and (3 = 0.76, each with 
isotropic angular distribution at the detector. 

Energy loss of the magnetic monopoles as they pass 
through the ice is modeled using the Bethe-Bloch formula 
as adapted by Ahlen [9] . The typical energy loss in ice is 6- 
10 GeV/cm. The light output and propagation is modeled 
by a version of PHOTONICS [10] specifically generated 
to work with Cherenkov cone angles associated with the 
different speeds simulated. The default light amplitude in 
the simulation is for a muon and this is scaled up using the 
formula of Tompkins [7]. 

The principal background to relativistic magnetic 
monopole searches with IceCube consists of down- 
going high energy atmospheric muon bundles induced by 
cosmic rays. Using a 2-component model, which assumes 
the cosmic ray flux to be composed of only protons and 
iron, cosmic ray primaries are simulated in the energy 
range from 10"* GeV to 10^^ GeV, where the energy spectra 
of the components have been fitted to the KASCADE 
data [11 J. The muons are generated with the air-shower 
simulation package CORSIKA [12] and handed to the 
detector simulation software. 

4 Search Strategies 

The general procedure used in the two following magnetic 
monopole searches is that of a blind analysis. Hence the op- 
timization of the data selection is based on simulated data, 
and only ^ 10% of the experimental data, referred to as 
the burn sample, is used for verification purposes. For the 
monopole signal, a flux of 5- 10~^^ cm~^ s~^ str~^ is used, 
roughly representing the lowest limits set by BAIKAL [13] 
and AMANDA [14]. 
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Figure 1 : Monopole flux limits set by the 22 string analysis 
versus beta, shown together with the previous best limits 
set by AMANDA [14] and BAIKAL [13]. The limits ap- 
ply to an isotropic flux at the detector. Also shown are the 
expected sensitivity of the 40 string analysis, limits set by 
flie MACRO [18] experiment and the Parker Bound [19]. 

4.1 IceCube-22 

The first search for relativistic magnetic monopoles was 
performed on the 2007 data run with the detector operat- 
ing with 22 strings representing a volume of ^ 0.3 km^. 
The selection strategy focused on the brightness of the 
monopole by first defining a hit with a much higher thresh- 
old than the actual DOM, set to be when the fADC wave- 
form reached a saturation level of 1022 ADC counts. This 
reduced hits from background to be ~ 10 m away from the 
particle track while signal hits were ^ 10-60 m away from 
the monopole track. Once the sequence of hits was defined, 
a simple analytic reconstruction was performed based on a 
plane wave hypothesis and minimizing the variance relative 
to the actual hit positions, termed 'linefit' . The result of us- 
ing the higher threshold for a hit leads to an improvement 
in angular accuracy from 5-10 degrees for all hits down to 
~ 2 degrees. 

The analysis aimed to enhance the sensitivity to slower 
monopoles by binning the data based on speed reconstruc- 
tion. UnbUnding revealed that the background simulation 
was not reproducing the tail of the speed distribution where 
the slower signal was expected to be. This resulted in ob- 
vious background events surviving into the final sample. 
Comparisons based on additional characteristics of wave- 
forms and spatial extent of the hits clearly distinguished 
the events as background muon bundles. After determining 
no monopole events were recorded, the following analy- 
sis was based only on improved simulation rather than the 
bum sample or other experimental data. The only changes 
involved a shght tightening of quality cuts motivated by the 
new simulation and abandoning any speed biiming. 

The final cut is in the plane of reconstructed zenith direc- 
tion and the number of bright hits. For up-going events. 
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Figure 2: Ratio of the number of hit DOMs and total 
number of photo-electrons for simulated atmospheric muon 
background (dashed line), 10% experimental data (black 
markers) and simulated monopole signal (grey histograms). 
A precut (Nch > 50) has been applied to reduce the discrep- 
ancy between experimental data and Monte-Carlo simula- 
tion. 

that are largely background free, a simple cut on the num- 
ber of hits is applied. For down-going directions, where 
muon bundles produced by high-energy cosmic rays domi- 
nate, the cut increases in strength as it approaches the verti- 
cal. The cut is optimized on background and signal Monte 
Carlo by considering the combination that minimizes the 
model rejection factor [15] for an isotropic and monoener- 
getic monopole flux. 

No data events survived the final cut and upper limits were 
placed on the flux of magnetic monopoles for speeds of 
(5 > 0.8, which represent a factor of 10 improvement over 
the AMANDA analysis [16]. Figure 1 shows the limits, as 
a function of beta, that result for an isotropic flux at the 
detector. For a complete description of the analysis as well 
as how these limits transform to an isotropic flux at the 
Earth's surface, see [17]. 

4.2 IceCube-40 

A relativistic magnetic monopole search is also currently 
underway using the data taken with the 40-string detector. 
The selection strategy remained focused on the brightness 
of the monopoles, however, a different observable was cho- 
sen. The primary parameter used was the ratio of the to- 
tal number of photo-electrons (NPE), which is estimated 
by unfolding the f ADC waveforms, and the number of hit 
DOMs (Nch). While Nch and NPE are in principle in- 
tegers, their values are large for bright events and there- 
fore the ratio provides a wider range of values than us- 
ing the number of saturated DOMs. For bright monopoles 
NPE/Nch is expected to be larger than for atmospheric 
muon bundles, since the number of DOMs within a ra- 
dius r around a particle track is proportional to r^, while 
the number of photons decreases as exp(— r/A), where A 



Figure 3: Final cut on time over threshold to NPE ratio. 
The cut is shown as the thick vertical line at 1.6. 

is the effective absorption length. Figure 2 shows the dis- 
tribution of NPE/Nch after a precut (Nch > 50) has been 
applied to remove events where a particle track outside the 
instrumented volume triggers the detector. This also sig- 
nificantly improves the agreement of experimental data and 
Monte-Carlo simulation in the signal region. 

The actual cut was made in the plane of reconstructed di- 
rection and the inverse of the NPE to Nch ratio, where di- 
rection is reconstructed with the linefit method using all hit 
DOMs. In analogy to the previous search we made a sim- 
ple cut on the Nch to NPE ratio for up-going events and 
a cut increasing in strength toward vertical down-going di- 
rection. The downside of choosing these cut parameters is 
that monopoles with (3 = 0.76, i.e. close to the Cherenkov 
threshold, cannot be effectively discriminated from back- 
ground. The reason for that is the strong decrease of emit- 
ted Cherenkov light as the speed approaches the threshold. 
This narrows the speed interval where the analysis is sen- 
sitive, however, only by a small fraction. As an example 
figure 4 shows the 2-dimensional distribution of the cut pa- 
rameters for monopoles with /3 = 0.8 and simulated back- 
ground. 

The parameter of the final cut is defined as the ratio of the 
accumulated time the fADC waveforms stay above a cer- 
tain threshold (ToT) and again the total number of photo- 
electrons. Figure 3 shows the distribution of the cut parme- 
ter. The cut value is set so that no data events survive the fi- 
nal cut. The estimated sensitivities to magnetic monopoles 
wifli speeds /3 = 0.995, = 0.9 and /3 = 0.8 at the de- 
tector are given in Table 1 and are also shown in Figure 1. 
Though no optimization has been applied to the cut values 
so far, the sensitivity is improved over the analysis using 
the 22 string detector by a factor of ~ 4. 

Looking to the future, visual inspection of burn-sample 
events has revealed two classes of background events 
which are evading the current cut conditions. The first class 
is events where a single DOM registers a signal much larger 
than the remaining DOMs. These events are beheved to 
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Figure 4: Distribution of Nch/NPE vs Cos(d) for monopoles with /3 = 0.8 (left) and for simulated muon bundles (right). 
A zenith angle of = ° corresponds to vertically down going. The cut is shown as a grey line, excluding the upper area. 
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Table 1: Effective area, expected signal and sensitivity 
(90% C.L.) at the detector for a full year of data. 



occur when a particle experiences a stochastic energy loss 
very close to a DOM. The second and more frequent event 
class appears to be atmospheric muons with high inclina- 
tion grazing the comers of the detector resulting in a hit 
pattern easily leading to misreconstruction. The occurrence 
of these events may be partially caused by the asymmetry 
of the 40 string detector. Even though no data events sur- 
vived the last cut, the final version of the analysis will likely 
include safety cuts for these event classes. 

5 Conclusion 

The analysis performed on the data taken with the 22 string 
IceCube detector has improved the limits on the flux of 
magnetic monopoles with speeds 13 > 0.8 by an order of 
magnitude. They are presently the most stringent experi- 
mental limits. The ongoing analysis of data taken with the 
40 string detector will likely improve these limits by an- 
other factor of ~ 4. However, these results are preliminary 
and wiU be refined. 
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